function process_convergence_2d(vector_error, vector_N)
error_var_name = inputname(1);
convergence_rate = zeros(size(vector_N,2), 1);
for i = 2:size(vector_N,2)
    convergence_rate(i) = log(vector_error(i-1)/vector_error(i))/log(vector_N(i)/vector_N(i-1));
end
T = table(compose("%5.0f", vector_N'), compose("%5.4e", vector_error), compose("%5.5f", convergence_rate),...
    'VariableNames', {'N', error_var_name, 'Rate'});
disp(T);
end